Dynamics of Vortex Formation in Merging Bose-Einstein Condensate Fragments 
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We study the formation of vortices in a Bose-Einstein condensate (BEG) that has been prepared 
by allowing isolated and independent condensed fragments to merge together. We focus on the 
experimental setup of Scherer et al. [Phys. Rev. Lett. 98, 110402 (2007)], where three BECs are 
created in a magnetic trap that is segmented into three regions by a repulsive optical potential; the 
BECs merge together as the optical potential is removed. First, we study the two-dimensional case, 
in particular we examine the effects of the relative phases of the different fragments and the removal 
rate of the optical potential on the vortex formation. We find that many vortices are created 
by instant removal of the optical potential regardless of relative phases, and that fewer vortices 
are created if the intensity of the optical potential is gradually ramped down and the condensed 
fragments gradually merge. In all cases, self-annihilation of vortices of opposite charge is observed. 
We also find that for sufficiently long barrier ramp times, the initial relative phases between the 
fragments leave a clear imprint on the resulting topological configuration. Finally, we study the 
three-dimensional system and the formation of vortex lines and vortex rings due to the merger of 
the BEC fragments; our results illustrate how the relevant vorticity is manifested for appropriate 
phase differences, as well as how it may be masked by the planar projections observed experimentally. 
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I. INTRODUCTION. 



The formation, stability and dynamics of vortex-like 
structures has been a long-standing theme of interest in 
many areas of physics, including classical fluid mechanics 
[l| , superfluidity and superconductivity 0, H, 3 , and cos- 
mology 0]. Moreover, in the past decade, there has been 
a tremendous growth of excitement in this topic in the 
branches of atomic and optical physics. This has been 
propelled by considerable experimental and theoretical 
advances in the fields of nonlinear optics and Bose- 
Einstein condensates (BECs) in dilute alkali vapors 0, Q 
(see also Ref. 

Focusing more specifically on the rapidly growing area 
of BECs [lOl , one can recognize that the study of vortices 
has been central to the relevant literature. In particular, 
as concerns the experimental efforts, the original obser- 
vations of single [ill, and multiple vortices [ll] was 
soon followed by the realization of robust lattices of large 
numbers of vortices [T3 |. Subsequent studies turned to 
higher-charged structures such as vortices of topological 
charge S = 2 and even 5" = 4 [l^ and illustrating dynam- 
ical instability of these topological objects [l^. On the 
other hand; the abundance of experimental results has 
stirred an intense theoretical interest in the conditions 
under which such vortices and vortex lattices would be 
robust and experimentally observable. Most often, vor- 
tex existence and stability issues were examined in the 
framework of the standard parabolic confining potential 
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(typically produced by magnetic traps). In that frame- 
work, and in the two-dimensional (2D) case, vortices of 
charge S = 1 were found to be stable, while vortices of 
higher charge {S = 2, 3) were shown to be potentially 
unstable [l7j (depending on the atom numbers). Later, 
similar results were found for vortices of S* = 4 p^ . 
while the availability of more substantial computational 
resources has more recently led to similar conclusions in 
the fully three-dimensional (3D) case [l^, [2^ . The stud- 
ies of Rcfs. [2lj and [S^l examined the various scenarios 
of break-up of higher-charge vortices during dynamical 
evolution simulations for repulsive and attractive inter- 
actions respectively. Furthermore, Ref. [2^ considered 
such vortices riding on the background of not just the 
ground state, but also of higher, ring-like, excited states 
of the system. It should also be mentioned that these 
advancements have motivated the development of math- 
ematically rigorous tools in order to study the spectrum 
of such vortex modes. Such methods involve the use of 
the Evans function [23|, or the use of the index theorem 
evaluating the number of potentially unstable eigendirec- 
tions Hi]. 

Although the existence and stability of fundamental 
and higher charge vortices has been examined extensively 
as indicated above, the formation of such vortex struc- 
tures is far less studied. In particular, while seminal in- 
terference experiments (demonstrating that BECs are co- 
herent matter waves) were reported as early as a decade 
ago [2^, the role of interference between BECs in vor- 
tex generation was experimentally studied only recently 
(27| . This work proposed and examined the formation 
of vortices resulting from the interference and controlled 
merging of three condensed fragments, where the frag- 
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ments were essentially independent BECs separated by 
an optical potential barrier. Such a process has close 
ties to elements of topological defect formation in phase 
transitions, as proposed by Kibble Q and Zurek In- 
terestingly, the experimental work was almost concurrent 
with a theoretical study investigating a sinipler elongated 
barrier separating two independent BECs [28| ; in the lat- 
ter setting, the interference forms a dark soliton whose 
bending and subsequent breakup due to the manifesta- 
tion of the transverse modulational instability also result 
in vortices. 

In the present work, we expand on these considerations 
and study in detail, by means of systematic numerical 
simulations, the formation of vortices in a setting closely 
matching the one of the experiment in Ref. [27[. Our 
purpose is to get a deeper insight into this interference- 
induced vortex formation mechanism, investigating fun- 
damental features, such as the number and lifetime of 
ensuing vortices. This is done upon studying in de- 
tail the parametric dependences influencing the rele- 
vant experimental observations. Specifically, we quan- 
tify the above mentioned features as functions of the 
elimination/ramping-down time of the optical barrier be- 
tween the fragments, or the initial relative phases be- 
tween the original independent fragments. Our investi- 
gation chiefly refers to a 2D setting, but we also illustrate 
how the results are generalized in the pertinent 3D case. 
Notice that our considerations are motivated not only by 
their direct bearing on the experiments of Ref. [27| , but 
also by their relevance to studies of spontaneous symme- 
try breaking during phase transitions 0, . 

Our presentation is structured as follows. In Sec. [TTl 
we briefly summarize the setup of our computational ex- 
periments. In Sec. mil we study the interference of three 
BEG fragments in the 2D and 3D setup; special attention 
is payed at the role of the relative phases between the dif- 
ferent condensates, and the ramp-down time of the laser 
sheet barrier responsible for separating the three frag- 
ments. In Sec. IIVI we briefly comment on the relation 
between numerical and experimental results. Finally, in 
Sec. |Vl we summarize our findings, present our conclu- 
sions, and discuss some possible extensions of this work. 



II. SETUP 

We consider a BEG at a temperature close to zero, 
where quantum or thermal fluctuations arc negligible 
(note that finite temperature effects arc briefly discussed 
at the end of Sec. IIIip . This system can accurately be 
described by a mean-field theoretical model, namely the 
Gross-Pitaevskii equation (GPE) p^ : 



dip 



_ V2 + y(r;t)+g|Vf 
2m 



(1) 



where ip = ip{r,t) is the condensate wavcfunction (with 
n{r,t) = |?/;(r,t)|^ being the atomic density of the con- 
densate), m is the atomic mass, the coupling constant 



g = ATrh?as/m measures the strength of inter-atomic in- 
teractions and as is the s-wave scattering length. The 
potential V{r; t) in the GPE is taken to be of the form, 



V{r;t) = VuT{r) + a{t)VL{r) 



(2) 



where the two components in the right-hand side of 
Eq. ([2]) are a harmonic magnetic trap, Vmt{^) ~ 
^m{ijj'^x^ + uj^y^ + ojIz^) , with trapping frequencies lOx = 
bjy = 211 X 7.4 Hz and = 27r x 14.1 Hz, and the three- 
armed time-dependent optical barrier, a(i)VL(r), used in 
the experiments of Ref. [23] . This three-armed potential 
induces a separation of the ground state of the condensate 
into three different fragments (see top-left and top-center 
panels in Fig. [T] and Fig. [TlT a)). Note that the function 
a{t) in Eq. ^ describes the ramping down of the optical 
barrier. The maximum initial barrier energy for the po- 
tential is taken to be ao = a(0) = 26 /cb nK [23|, where 
/cs is Boltzmann's constant. 

In the numerical simulations, the ground state of the 
system is obtained by relaxation (imaginary time inte- 
gration) initiated with the Thomas-Fermi (TF) approxi- 
mation (see top-center panel in Fig. [T]) 



V'(r, 0) = Vmax{0,At- V"(r)} x (^(r). 



(3) 



where /i = Sfcs nK [23| is the chemical potential, and 
(^(r) contains the chosen phase of the different wells sep- 
arated into the three regions A, B and G depicted in the 
top-right panel of Fig. [1] Let us denote by </>! , 4>2, and 
(/)3 the initial phases in regions A, B and G, respectively. 
It is important to note that while in the experiments of 
Ref. [27 1 the different fragments have uncorrelated (ran- 
dom) phases, in our numerical simulations we arc able to 
control their initial phases and, more importantly, their 
relative phases. Our numerical experiments, emulating 
the experimental sequence of Ref. j27| . are performed 
with the output of imaginary time relaxation used as 
initial condition for the full dynamics of Eq. (HI). The 
bottom row of Fig. [T] depicts the evolution of the phase 
for the case (pk ~ 27rfc/3 during relaxation. As can be 
observed from the figure, the initial "seed" {t = 0) starts 
with sharp phase boundaries around the localized frag- 
ments. As the relaxation procedure evolves, the phase 
boundaries become smoother and adjust to the bound- 
aries for the top-right panel of the figure. The phase 
profile seems to settle after 10 ms of imaginary time re- 
laxation. In our simulations, 50 ms of imaginary time 
relaxation was used to ensure proper convergence to the 
steady-state solution. 



III. NUMERICS. 
A. Two-dimensional BECs. 

For the 2D rendering of the experiment of Ref. [27| . 
we restrict our system to the (x,y) coordinates and we 
use the same chemical potential as in Ref. [23|. First we 
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FIG. 1: (Color online) Top-left: intensity profile of the op- 
tical potential responsible for segmenting the potential well 
into three local minima. Top-center: Thomas-Fermi approxi- 
mation used as an initial condition for our relaxation method 
to obtain the ground state of the system. Top-right: regions 
A, B, and C with respective phases (jti, 4^2, and (^3. The bot- 
tom row of panels depicts the evolution of the phase during 
our relaxation (imaginary time relaxation) towards the initial 
steady state with different phases for the fragments. This ex- 
ample shows the phase for the case (j)k = 27rfc/3 at the times 
indicated. In all panels, the field of view is approximately 70 
/im per side. The axis numbers indicate x and y coordinates 
relative to the center of the unsegmented harmonic trap. 

explore the effect of the ramp-down time of the potential 
barrier on the formation of vortices through the merging 
of the different fragments of the condensate. For this 
purpose we use a linear ramp: 

a(i) -max 1^(4- i), o|, (4) 

where ao = a(0) is the maximum barrier energy as de- 
fined above and tb is the ramping time (in ms) of the 
barrier; note that a similar ramp was used in the experi- 
ments of Ref. [13 . 

We monitor the formation of vortices, as well as the 
overall vorticity of the system, using various diagnostics. 
These are based on the corresponding fluid velocity of 
the superfluid given by [s^l , 

ih ip*Vip - ipVip* 
v^- = -7; r-72 > (5) 

where (•)* stands for complex conjugation. The fluid vor- 
ticity is then dcflncd as a; = V x v,,. The results for the 
merger of the three BEC fragments with relative phases 
(j)k = 27rfc/3, k e {1, 2, 3}, for different merging times are 
depicted in Fig. O As can be seen from the figure, the 
number of vortex pairs nucleated by the merger is ex- 
tremely sensitive to the ramping time tb- Shorter ramp- 
ing times give rise to an extremely rich vorticity pat- 
tern as the fragments merge [see, for example, the top 
two rows in Fig. [2] corresponding to the instantaneous 




FIG. 2: (Color online) Evolution of the 2D condensate den- 
sity (respective top series of panels) and vorticity (respective 
bottom series of panels) for three different ramp-down times 
of the optical potential barriers. From top to bottom, the 
three sequences correspond to tb = ms (first and second 
row), ti, = 25 ms (third and fourth row), and tt = 50 ms 
(fifth and sixth row). The times are indicated in the panels 
in ms and the field of view is approximately 70 ^m per side. 
For the initial conditions in all cases, the different condensed 
fragments have relative phases of 27r/3, namely, (jtk ~ 27rfc/3 

(fee {1,2,3}). 



(tb = 0) removal of the barrier], including the appear- 
ance of structures resembling vortex streets. However, 
for longer ramp times, i.e. slower ramping, just a handful 
of vortices are nucleated. In fact, for tb > 100 ms (re- 
sults not shown here), the only vortex that is nucleated is 
the central one. It is evident that, independently of the 
ramping time, the central vortex is always formed for this 
pair of relative phases between the fragments. This vor- 
tex is the consequence of the intrinsic vorticity present in 
the initial condition where the three fragments have been 
phase imprinted with a total of a 27r phase gain about the 
condensate center. 

In Fig. [3] we present a spatio-temporal rendering of the 
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FIG. 3: (Color online) Evolution of the vortex structures in 
the 2D condensate density for tt = 25 ms and (f>k = 2nk/3 
(cf. middle rows in Fig. [2J . 



vortex formation for the middle row example of Fig. [2] 
(i.e., phases given by (f>k = 27rfc/3 and a ramp down 
time of tb = 25 ms). In the figure we depict a space- 
time contour plot of the vorticity where blue/red con- 
tours correspond to negatively/positively charged vor- 
tices, respectively. The figure clearly shows the forma- 
tion of pairs of vortices with opposite charge, some of 
which self-annihilate at later times, while others oscillate 
together with the cloud: upon formation, they expand to 
the rims of the cloud and are subsequently reflected from 
its outskirts and contract anew, following the density pro- 
file oscillations resulting from the merging process. 

In order to measure the vorticity generated during the 
merger of the condensate fragments at any given time, 
we compute the expectation value of the z-component of 
the angular momentum of the BEC 



= ih / V* dei) dr = {tp, (r x p)zi') ■ 



(6) 



where f and p denote the corresponding position and 
momentum operators, the subscript z indicates the com- 
ponent of the corresponding cross-product and (•, •) indi- 
cates the complex inner product defined from CxC C. 
Using Eq. ^ and Eq. ([T|), we can evaluate the time 
derivative of this expectation value: 



dt 



ait) 



\i'Wrxp),VLir)dr. 



(7) 



In obtaining this result, we have assumed the isotropy 
of Vmt{i^) in the (x,y) plane. This result also has some 
important consequences including the fact that during 
the ramp down of the optical barrier, we should not ex- 
pect the angular momentum of the BEC to be conserved, 
while we should expect such a conservation to appear 
once a(t) = or, more generally, when the full potential 
is azimuthally isotropic in the {x,y) plane (cf. Figs. [3] 
and [6]). In the left column of Fig. [4] we depict the z- 
component of the angular momentum. The middle col- 
umn of the same figure shows a quantity that we refer to 
as total fluid velocity, defined as: 
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FIG. 4; (Color online) Vorticity indicators for the cases pre- 
sented in Fig. [2] (namely, tj, = ms, ti, = 25 ms, and ti, — 50 
ms, from top to bottom) with (j)k = 27rfc/3. The left panels 
correspond to the total angular momentum [Eq. ([S}] normal- 
ized by h, while the middle and right panels correspond to 
the total fluid velocity [Eq. ((Sjl] for the whole cloud (middle) 
and the central portion (see text) of the cloud (right). 



where V = / dr is the total volume of integration. Fi- 
nally, the right column of Fig. 3] shows the same quantity, 
as defined in Eq. ([5]), but only for the central portion of 
the cloud. The central portion of the cloud was defined 
as a square, centered at the center of the magnetic trap, 
with a side equal to 10% of the integration domain. This 
area corresponds approximately to the void area between 
the three initial fragments (see top-center panel in Fig.[T]). 

The three diagnostics defined above are shown in Fig. 2] 
for the same cases presented in Fig. [51 It is interesting 
to note how the presence of a time-dependent compo- 
nent in the potential generates angular momentum in 
the system. As can be observed in the figure, the total 
angular momentum increases (in absolute value) through 
the duration of the barrier ramp-down and then settles 
to a constant value that is larger (in absolute value) for 
slower ramps. Also, it is worth noting that the initial 
value of the angular momentum is different from zero 
[Lz{t = 0) ~ —800] due to the intrinsic vorticity car- 
ried by the out-of-phase fragments (see also discussion 
below). For our initial condition, most of the total an- 
gular momentum is seeded in the weak overlapping re- 
gion between the fragments where the phase gradient is 
large. What we can observe about the second and third 
diagnostics by comparing (the second row of) Fig. [2] and 
Fig. [3] with (the second row of) Fig. [3] is roughly the fol- 
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FIG. 5: (Color online) Same as in Fig. [2] for three different 
relative phases and a ramp down of tt = 25 ms. The top two, 
middle two and bottom two series correspond, respectively, 
to: a) <j)k = 0, b) = 0, 02 = 27r/3, and (j>3 = 47r/3, and c) 
(j>i = 0, (j>2 = 7r/3, and 4>3 = 27r/3. 



lowing: the integrated (throughout the cloud) velocity 
appears to peak when the filamentation in the pattern of 
Fig. [3] is maximal, e.g., around times of 50 and 90 ms. 
On the other hand, the same diagnostic integrated within 
the central core of the cloud peaks substantially earlier 
when the vortices are formed through the collision of the 
fragments around the end time of the ramp (i.e., around 
25 ms). Subsequently, and as the vortices are advected 
away from the core, the latter quantity decreases. 

We now turn to the examination of the effect of the 
relative initial phases between the different fragments. 
As an example, we show in Fig. [5] the behavior of the 
cloud density and its vorticity for three different phase 
combinations for a fixed ramp-down time of tb = 25 ms. 
As can be noticed, the complexity of the vorticity field is 
similar for the three cases shown in the figure. However, 
the only phase combination of the ones shown here that 
produces a vortex at the center of the cloud corresponds 
to 0fc = 27rfc/3, which can be understood by the intrin- 
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FIG. 6: (Color online) Vorticity indicators for the cases de- 
picted in Fig. [S] (the different panels are presented in the same 
manner as in Fig. 3)). 



sic vorticity already present in the initial configuration 
(see explanation below) . In Fig. [6] we depict the vortic- 
ity indicators for the three phase combinations of Fig. [S] 
As can be seen in the figure, the angular momentum Lz 
for 0fc = (top-left panel) and (j)k = nk/S (bottom-left 
panel) suffers almost no change since the initial config- 
uration does not carry any intrinsic vorticity and thus, 
its interaction with the ramping down barrier does not 
produce angular momentum. However, for = 27rfc/3 
(middle-left panel), as explained before, the system does 
gain angular momentum during the barrier removal. The 
total fluid velocity indicators (middle and right columns 
in Fig. [5] behave similarly for the three phase combina- 
tion with the notable difference that for the = 27rfc/3 
(middle-left panel) case, its final value is different from 
zero since the fragments produce a vortex at the center 
of the trap (due to the intrinsic vorticity carried by the 
initial condition). 

In order to follow in more detail the formation of vor- 
tices in the central portion of the cloud as a function of 
the relative initial phases of the different fragments, we 
perform systematic simulations for a large set of relative 
initial phases. The results are presented in Fig. [71 where 
we show the existence of vortices as a function of (02, '/'s) 
for = 0. Darker shades correspond to the presence of 
more vortices. The top panel of Fig. [7] corresponds to 
an immediate release of the barrier (ramp down time of 
tb = ms) where the presence of 0, 1, 2 or 3 vortices 
(white, light gray, gray, and black respectively) can be 
observed for different phase combinations. The middle 
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pa. l32| for reviews). In particular, each fragment can be 
thought of a "unit" and the whole configuration corre- 
sponds to a discrete vortex with three units with a net 
vorticity different from zero. The middle panel of the fig- 
ure clearly reveals that only within an arc of the second 
(and essentially symmetrically of the fourth) quadrant of 
the plane of the phases (02 ~ tt, 03 — tt), a single vortex 
will form at the center of the configuration. It is inter- 
esting that discrete vortex-like configurations consisting 
of three fragments have been considered in the context 
of nonlinear optics in Ref. (ssj and even in genuinely dis- 
crete systems as, e.g., in Ref. [13]. In the bottom panel 
of Fig. [7] we depict the difference between the top and 
middle panels to show the amount of vortices that are 
created exclusively by the fragment collision and not by 
the intrinsic vorticity of the initial configuration. It is im- 
portant to mention that some of the vortices counted in 
the top panel come from vortices that are created outside 
the central region but that migrate towards the center as 
time progresses. 



B. Three-dimensional BECs. 
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FIG. 7: Phase diagram for vortex formation at the trap center 
region as a function of the relative phases of the merging 
fragments (where 01 is assigned the value of 0). Top: Ramp 
down time of tt = ms (i.e., instantaneous removal) of the 
barrier. The different shades of gray indicate the amount of 
vorticity contained in the central region after 100 ms. White, 
light gray, gray, black correspond, respectively, to 0, 1, 2, and 
3 vortices. Middle: Same as above for a ramp down time of 
tb = 50 ms. This case only produces no vortices (white) or one 
vortex (gray). Bottom: difference between top and middle 
diagrams. This corresponds to the vortices formed by the 
collision of the fragments and not by the intrinsic vorticity of 
the initial configuration (which depends on the relative phases 
of the different fragments). 



panel depicts the same diagram, but for a ramping-down 
time of tb = 50 ms. It is clear that, for this relatively 
slow ramp down, the formation of a vortex in the cen- 
tral region is exclusively determined by the vorticity of 
the initial configuration. Namely, if any of the relative 
phases is larger than tt, the initial configuration resem- 
bles more that of a discrete vortex (see Refs. (sij and 



As shown in the previous section, the 2D setting lends 
itself to a more detailed examination of important fea- 
tures such as the parametric dependence on the ramping 
time and the relative phases of the fragments. Neverthe- 
less, it is important to also consider some of the delicate 
points particular to the 3D nature of the experiments 
and the observable quantities available within the exper- 
imental images. For this reason, we now focus on the 
3D setting, presenting results of the simulations relevant 
to the experiment of Ref. The setup is the same 

as in the previous section but we now use the full 3D 
space with the same chemical potential as before. Typ- 
ical results are shown in Figs. [5)1101 for a ramping down 
time of the potential barrier of 25 ms. Figure [5] corre- 
sponds to the case of equal initial phases 0^ = 0, while 
Fig- El pertains to initial phases 0^ — iirk/Z, and Fig. [TUl 
to phases (pk ~ 7rfc/3. The figures depict contour plots 
of the density (top rows) and vorticity (middle rows), 
as well as a z-projection of the density (bottom rows) 
as it would be observed in the laboratory. The vortex 
structure is considerably more complex in the 3D sce- 
nario because the vorticity does not show up as straight 
vortex lines but rather as a complex web of vortex fila- 
ments in various directions. As in the 2D case, there is 
the formation of a vertical vortex line at the center of 
the cloud (cf. second row in Fig. for the appropriate 
relative phases of the different fragments of the conden- 
sate with the same conditions as before. Nonetheless, in 
the 3D case, the central vortex line is prone to bending 
as can be clearly seen in the later stages of the dynam- 
ical evolution presented in the second row of Fig. [9l In 
fact, the vortex bending is even clearly visible in the z- 
projection (see third row of Fig.lH]). These 3D numerical 
experiments arc quite revealing in that the laboratory ex- 
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t=15 ms 
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FIG. 8: (Color online) Vortex formation by the merging of three-dimensional BECs. Top row: contour surfaces of constant 
atomic density. Middle row: contour surfaces of the corresponding absolute value of the vorticity. Bottom row: z-projection 
of the density distribution (i.e., column density along z) as it would be observed in the laboratory. The snapshots are taken 
at the indicated times (in ms) for an initial phase distribution corresponding to = and a ramping down time of 25 ms. 
For the contour-surface images, the axis on the left side of each plot represents the vertical (z) direction, expressed in units of 
microns. 



periments can only show projections of the density and 
thus missing to a great extent are the intricate vortex line 
dynamics. Importantly also, between the bending effect 
and the integrated view used in the experimental images, 
it is possible for the presence of vortcx-likc structures or 
filaments to be blurred (as in the later stages of Fig. [5]) or 
entirely lost (as in the later stages of Fig. [5] and Fig. [TU|) . 

As can be seen from Figs.lSHTUl the vorticity emerges at 
the early stages of the merger (t < 20 ms), through vor- 
ticity sheets that nucleate some of the vortex line struc- 
tures. Nonetheless, it is interesting to note that most of 
the vorticity is carried by vortex lines and vortex rings 
that are horizontal (except the notable case of the ver- 
tical vortex line depicted in the second row of Fig. [5]). 
This fact is also quite visible in the density contour plots 
for t = 60 ms where the horizontal vortex lines "pinch" 
the cloud and create peripheral horizontal ridges around 
the cloud. It is also possible to observe some vorticity in 
the bulk of the cloud that docs not directly come from 
the phase differences between the initial fragments, but 
from the actual turbulence that is created by the frag- 
ment collision. As an example, two small vortex rings 
are clearly visible in Fig. [5] for t — 60 ms (one close to 
the top of the cloud and the other one 1/3 from the bot- 
tom). We would like to stress the difficulty of capturing 
the vorticity at the edge of the cloud (where most vortic- 
ity is actually observed) in our numerical experiments. 
This is due to the fact that the vorticity is defined as 



the curl of the fluid velocity of Eq. ^ that is normalized 
by the density. The numerical effect is that close to the 
periphery (where the density is small) the fluid velocity 
corresponds to the ratio of small numbers which imposes 
great numerical difficulties. Nonetheless, by using a fine 
grid of 301 x 301 x 121 we are able to capture most of 
the delicate vorticity dynamics at the periphery of the 
cloud consisting, mostly, of horizontal vortex lines that 
are parallel to the periphery of the cloud. 

Another interesting phenomenon is the oscillation of 
the atomic cloud. The cloud starts with a larger horizon- 
tal extent compared to the vertical one and after merger 
creates an almost spherical cloud, which in turn elongates 
again in the horizontal direction after the fragments go 
"through" each other. This behavior repeats a few times 
until the cloud takes an approximate spherical shape (re- 
sults not shown here). 

We have also monitored the effects of damping due to 
the coupling of the condensed atoms to the thermal cloud. 
Equation ([1]) is obtained by supposing a dilute Bosc gas 
at a temperature close to absolute zero. However, at finite 
temperatures, but still smaller than the critical tempera- 
ture Tc for condensation, a fraction of the atoms are not 
condensed and form the so-called thermal cloud. In turn, 
this thermal cloud induces a damping on the dynamics 
of the condensed cloud. We used the approach of phe- 
nomenological damping [s^ described in Refs. [H, [s^l 
that relics on replacing the i in front of the time deriva- 
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FIG. 9: (Color online) Same as in Fig. [H]for an initial phase distribution corresponding to (f)k = k27v/S. 
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FIG. 10: (Color online) Same as in Fig. [8]for an initial phase distribution corresponding to (pk = fcvr/S. 



live in Eq. ([T]) by (i— 7), where 7 is the damping rate, and 
by renormahzing the solution at each iteration to keep 
the initial mass (number of atoms) constant during inte- 
gration. We tested values of 7 in the interval [0.01,0.1] 
that contains the value of 0.03 estimated in Ref. [Ill for a 
temperature T — O.lTc. The results of the phenomeno- 
logical damping are, qualitatively, very similar (results 



not shown here) to the effects of ramping down the bar- 
rier over longer time scales: larger damping resulting in 
a stronger suppression of the vorticity generated by the 
merger of the different cloud fragments. 



9 



(a) ^ 






(d) 






• 


# 


50 fim 


100 ftm 


100 urn 


100 urn 



FIG. 11: (a) In-situ phase-contrast image of three BECs 
trapped in a triple-well potential, (b)-(d) Absorption im- 
ages of BECs after 56 ms of ballistic expansion. Each BEC 
was created by merging three BECs, as described in the text. 



IV. COMPARISON OF NUMERICAL AND 
EXPERIMENTAL RESULTS 

The simulations described in the present work were di- 
rectly aimed at developing a more thorough understand- 
ing of both the experimental results of Ref. [23| and the 
dynamics of vortex formation during BEC merging and 
collisions. In this section, we thus briefly discuss the 
observed similarities and differences between the experi- 
mental and theoretical results. 

In the laboratory experiment, an optical potential was 
used to segment a harmonic trap before condensation 
was achieved; with additional evaporative cooling, three 
initially isolated and mutually independent (i.e., uncor- 
related phases between the different fragments) conden- 
sates were created. A phase-contrast image of three such 
BECs is given in Fig. [TlT a). which can be directly com- 
pared with the simulated data of Fig. [8] (bottom row, left 
image). The initial condition of incoherent condensate 
fragments serves as the conceptual basis behind the mo- 
tivation to impose and examine various relative phases 
between the fragments in this work's simulations. In the 
experiment, the optical barrier was ramped off approx- 
imately linearly over time scales between 50 ms and 3 
s, significantly longer than the time scales considered in 
the simulations. At the end of the merging process, im- 
mediately after the barrier was completely removed, the 
fully merged BEC was released from the harmonic trap 
and allowed to ballistically expand for 56 ms to enable 
imaging of vortex cores. Example images of merged and 
expanded BECs are shown in Fig. [TIT b) through (d). 
These images can be compared with the simulated data 
of Figs. [8l through [TOl (bottom rows); note that the sim- 
ulations do not involve an expansion stage. 

In both experimental results and in the simulations 
shown here, we note the following important similarities. 
First, it is clear that vortex cores may be formed during 
the merging process. Second, as noted in the experimen- 
tal work, the vortex formation process should depend 
upon relative phases between the condensates. This con- 
clusion is borne out by the present work. Specifically, 
the simple analysis of a slow merging process leading to 
a 25% probability of vortex formation, as presented in 
the experimental work, matches the results of the simu- 
lations summarized in Fig. [T] Finally, experimental and 
simulated results show that faster merging leads to more 



vortices initially created, but these vortices may self- 
annihilate with time by holding the fully merged BEC 
in the trap. 

However, there is one notable quantitative difference 
between the experimental and numerical results. In the 
experimental work, for barrier ramp-down times longer 
than 1 s, single vortices were experimentally observed 
in approximately 25% of the images obtained directly 
at the conclusion of the merging process. Multiple vor- 
tex cores were not observed under these conditions. For 
faster ramps, experimental images containing either sin- 
gle or multiple cores were more often obtained, with sig- 
nificantly more than 25% of the images containing at 
least one vortex core. The results of the numerical data 
show that the slow merging limit (where multiple vortices 
cease to be created during merging) is reached for merg- 
ing times that are much shorter than in the experiment. 
In other words, images with multiple vortices are seen 
in the experiment under conditions where the numerical 
results would suggest that a given BEC should have at 
most one vortex. 

There are a few possible sources of this discrepancy. At 
first glance, it might appear that the spontaneous forma- 
tion of vortices in BECs during evaporative cooling in 
an axi-symmetric harmonic trap, as noted in Ref. [27|, 
could play a role in higher percentage of vortices seen in 
the experiment. Such vortex formation processes can not 
be described by the GPE and are thus not observable in 
the simulations of this work. However, due to angular 
momentum damping and self-annihilation of vortices in 
the asymmetric local potentials in which the three BEC 
fragments grow, we believe that vortices that might be 
spontaneously created in one (or more) of the three BECs 
are unlikely to survive at rates that would affect the ex- 
perimental observations of Ref. [23| . This possible source 
for the quantitative discrepancy could be tested, for ex- 
ample, with simulation methods based on the stochastic 
GPE li. 

Perhaps a more likely source of the quantitative dis- 
crepancy may lie in the optical potential energy or shape; 
differences between the experiment and simulations re- 
garding barrier heights, widths, and ramp-down trajec- 
tory might induce more vortices to form during merging. 
For example, if center-of-mass oscillations of the cloud 
were induced in the experiment, atomic fiuid flow around 
the central portion of the optical barrier could induce 
formation of vortices and lead to increased vortex ob- 
servation rates in the experiment. Such processes could 
be studied in future GPE simulations in order to further 
characterize dynamical processes that may be involved 
in vortex formation. It might also be possible that im- 
perfections in the true optical barrier used in the experi- 
ment could pin vortices for a portion of the barrier ramp 
process, and significantly alter the vortex formation and 
annihilation process. Finally, much of the present work 
focused on the central portion of the trap, a region that 
encompasses the center of the merged BEC. In the ex- 
perimental work, however, vortices were most often seen 
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further from the BEC center. 



V. CONCLUSIONS 

We have studied the formation and subsequent evohi- 
tion of vortex structures and filaments in a system di- 
rectly simulating the experimental setup of Ref . [l^] ■ In 
particular, we have considered the case of three indepen- 
dent fragments (of variable initial relative phases) and 
how these merge upon the ramping down and eventual 
removal of the optical barrier that separates them. While 
there are many similarities between the numerical results 
and the experimental results of Ref. [131 , the numerical 
simulations importantly show features and new dynamics 
not discussed or observed in the experimental work. 

The first part of our study concerns the simpler two- 
dimensional setting, where it is straightforward to ob- 
serve the interference of the independent matter waves, 
the ensuing formation of vortices, as well as their mo- 
tion within the cloud, as a function of different param- 
eters such as the ramping-down time or the initial rela- 
tive phases between the fragments. Different diagnostics 
for the vorticity were developed in the process (such as 
the z-component of the angular momentum, or the inte- 
grated velocity of the flow throughout the cloud or near 
its center) and their dynamics was explained based on 
the evolution simulations. Principal findings of this part 
of the work included the formation of smaller numbers of 
vortices as the ramping-down time was increased and the 
formation of a single vortex in the core of the condensate 
for appropriate, discrete- vortex-like relations between the 
phases of the different fragments. 

The second part of our work explored how the features 
found in the two-dimensional setup are generalized in a 
fully three-dimensional setting, and how these affect the 
measurement process through, e.g., the projection of the 
BEC density on the plane. Key features of the latter 
dynamical evolutions involved the blurring of the vor- 
tex dynamics by the projection process coupled with the 
spontaneous vortex bending even when the different frag- 
ments have the appropriate phase relation to generate a 
vortex through their merging. In the 3D setting, the 
vorticity emerged in the form of vorticity sheets inducing 
vortex filaments (most often in a horizontal form) which 
led to pinching effects at the vortex cloud periphery and 
the formation of corresponding ridges in the atomic den- 
sity profile. 

Our work is related to a physical mechanism that has 
previously been discussed in the context of topological 
defect formation and trapping during phase transitions. 



often referred to as the Kibble-Zurek (KZ) mechanism 
d, Q • For the case of cooling of an atomic gas through 
the BEC phase transition [40|, the KZ scenario involves 
the growth and subsequent merging of phase-incoherent 
regions of the atomic cloud, with vortices being trapped 
in the merging process as the BEC grows in size and atom 
number. The full KZ mechanism involves physics beyond 
the scope of the simulations of this work; however, our 
simulations show dynamical processes that may relate to 
the portion of the KZ mechanism involving the merging 
of phase-incoherent regions of condensed atoms. Our nu- 
merical results showing the phase relationships involved 
in vortex formation during merging are consistent with 
basic notions of this portion of the KZ mechanism. 

There are many interesting questions to consider for 
future work in the present framework. Firstly, it is clear 
that further analysis of experimental data and further 
variations of simulation parameters will be needed in or- 
der to resolve the quantitative differences between exper- 
imental and numerical results, as discussed in See. IIVI 
Perhaps additional light on this question (and an inter- 
esting diagnostic in its own right) would be the exam- 
ination of the integrated density along the (a;, z) plane 
which should perhaps detect some of the horizontal vor- 
ticity filaments illustrated herein. A particularly chal- 
lenging (and more general) question along the same vein 
concerns the extent to which it may be possible to re- 
construct the fully three-dimensional cloud density from 
such projections. 

Both from an experimental and from a theoretical 
point of view it would be interesting to extend the present 
considerations also to multi-component condensates in 
order to examine the potential formation of vortex-like 
filaments and structures in settings similar to the ones 
presented herein (e.g., containing fragments from differ- 
ent components). In the latter setting, there would ex- 
ist an exciting interplay between the interference mecha- 
nisms and the formation of the coherent structures, and 
the phase separation dynamics between the components; 
see e.g., the recent experimental results of Ref. [32| and 
references therein. 
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